Bilayer graphene dual-gate nanodevice: An ab initio simulation 
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We study the electronic transport properties of a dual-gated bilayer graphene nanodevice via 
first principles calculations. We investigate the electric current as a function of gate length and 
temperature. Under the action of an external electrical field we show that even for gate lengths up 
100 A, a non zero current is exhibited. The results can be explained by the presence of a tunneling 
regime due the remanescent states in the gap. We also discuss the conditions to reach the charge 
neutrality point in a system free of defects and extrinsic carrier doping. 
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There is a consensus in the scientific community 
that the scaling down of silicon-based metal-oxide- 
semiconductor field-effect transistors (MOSFETs) is ap- 
proaching its limits. Therefore, there is a pursuit for dif- 
ferent materials to replace the silicon paradigm^ . Among 
the potential candidates, graphene has attracted the at- 
tention of research groups after the seminal work of the 
Manchester group^. Graphene is a truly two-dimensional 
(2D) material, with a honeycomb structure, and is a 
zero band-gap semiconductor. The valence and conduc- 
tion bands close to the Fermi level are cone shaped, 
with a linear energy-momentum relation at the two k 
points {K and K^) of the Brillouin zone, called Dirac 
points^. As a result, graphene has several properties that 
make it very interesting to use in the manufacture of de- 
vices, such as its very high electronic mobilities, up to 
200 000 cm^V~^5~^^'^, with electrons and holes behav- 
ing as massless fermions near the K and points. Thus, 
charge carriers can travel for micrometers at room tem- 
perature without any scattering. 

However, a drawback in using graphene as an elec- 
tronic material, such as in logic devices, is its lack of an 
energy gap. A number of different approaches have been 
proposed to open such a bandgap on graphene^^^. An- 
other way to create the gap is to use bilayer graphene, 
breaking the inversion symmetry through an external 
perturbation, such as the application of an electrical field 
perpendicular to the layers. Moreover, this band gap can 
be tuned by varying the field strength. This fact has been 
confirmed by both photoemission and optical absorption 
experimentsii — . A few theoretical studies have also in- 
vestigated this approach to open a band gap in bilayer 
graphene via the application of an electric field along the 
whole systemMii^. Recently, Avouris' group observed a 
transport band gap in biased bilayer graphene^^, also 
demonstrating its potential for applications in digital 
electronics. Thus, understanding the intrinsic properties 
of devices based on bilayer graphene in nano-dimensions 
is important for future developments in the area. 

In this Brief Report we present investigations of the 
electronic transport properties of a dual-gated bilayer 
graphene nanodevice, as a function of the gate length 
(Lgate) and temperature, via first-principles calculations. 




FIG. 1: Schematic view of a dual-gated bilayer graphene. 



We show an exponential decrease of the current as a func- 
tion of the gate length, which is characteristic of a tunnel- 
ing regime. As a result, the system does not reach a zero 
current. Also, we discuss the conditions to achieve the 
charge neutrality point (CNF) in a system that is not af- 
fected by defects and extrinsic carrier doping. We present 
results for electric fields comparable to the experimental 
electric displacement fields, l.OV/ nm < |£^| < 5.0V/ nm, 
for room temperature and A.bK. We use the Landauer- 
Biittiker model, with a Hamiltonian generated through 
ab initio density functional theory (DFT),^^ coupled with 
a non equilibrium Green's function formalism in a fully 
self-consistent procedure (NEGF-DFT)^^. To investigate 
the effect of finite gates we implement a real-space Pois- 
son solver with multigrid techniques'^. 

In Fig. [1] we show a schematic representation of the 
bilayer dual-gated graphene nanotransistor considered in 
our calculations. The atomic geometry is composed of 
an AB stacked bilayer graphene (BG) with a scatter- 
ing region up to 15 nm in length, sandwiched by left and 
right leads. A source-drain voltage {Vds) is applied across 
the system defining the current direction. The charge 
transport properties were calculated with the NEGF- 
DFT formalism, which uses the Kohn-Sham Hamiltonian 
and overlap matrices computed with the siesta code^^« 
to self-consistently obtain the density matrix. In order to 
investigate a finite dual-gate system, we have done two 
main modifications in the transampa code^': (i) The 
inclusion of k points that are transversal to the trans- 
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port direction {k±)^ where for each kj_ point in reciprocal 
space we find a /c^ -dependent density D{k±_)^y^ by means 
of an integration of the lesser Green's function^*^^ , and 
transmittance T{E^ k±)^ given by>^ 



Tr 
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G^^^\E^ k±) is the retarded (advanced) Green's function, 
whereas, Tj^-^j^^{E) are the terms that couple the leads to 
the scattering region. The total density and transmit- 
tance are obtained upon an integration in the reciprocal 
space with the following equations: 
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These modifications allow the simulations of truly 2D 
systems. Once the density matrix convergence is 
achieved, the current is calculated within the Landauer- 
Biittiker^^ model in the non-interacting approach of 
Meir-Wingreen^^, where the current is given by: 
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TiE)[f{E-^L)-fiE-m)]dE, (4) 



T{E) is given by ([3]), and f{E — jj^l/r) ^ire the Fermi- 
Dirac distributions for the left (L) /right (R) leads. We 
used the local density approximation (LDA) for the 
exchange-correlation functional^^ since it correctly de- 
scribes the graphene interlayer distance without the in- 
clusion of van der Waals corrections (ii) The sec- 
ond modification was the inclusion of a real-space Pois- 
son solver that allows non-periodic solutions for the 
Hartree potential, necessary to set different values for 
the top {Vtg) and back (Vbg) gate voltages with finite 
gate lengths. Therefore, we solve the Poisson equation in 
a rectangular box, fixing the upper and lower boundary 
conditions, where the values of Vtg and Vbg are defined 
only in a finite region of length Lgate (see Fig. [1]). In our 
investigation we consider Lgate = 1, 2, 3, 4, 5, and 10 nm, 
Lscatt ^ 15 nm, and a distance between the upper and 
lower boundaries of 20 A. 

We guarantee the convergence of the transmittance us- 
ing a double- C basis, 300 Ry for the mesh cut-off and 
800 k^ points in the Brillouin zone with the Monkhorst- 
Pack scheme^^. We used a fully relaxed geometry with 
a carbon-carbon distance of 1.428 A and an interlayer 
distance of 3.203 A. In the transport calculations, a tiny 
imaginary value of 10~^ Ry is added to the energy in 
order to calculate the retarded Green's function, and the 
electron density was calculated with 60 energy points and 
5 poles to evaluate the integral of the lesser Green's func- 
tion. 

The fiowing current Ids as a function of the source- 
drain voltage Vds is depicted in Figs. [2fa) and 2(b) for 




-i — ' — I — ' 

=+2.5V; V 



\.5N 



•-• No Gate 



300 
250 



200 
150^ 



100 



'50 



V,JmV] 



10 20 30 40 
V,JmV] 



FIG. 2: (Color online) In (a) we present Ids x Vds curves 
varying AV, for 300ii^ and Lgate = 5.0nm. In (b) Ids x Vds 
curves varying Lgate, for 300i^ and AV = 5.0V. 



room temperature. In our simulations the gate voltages 
Vbg and Vtg are tuned independently. In Fig. [2] all curves 
are for gate voltages with Vbg = — "^t^, where the volt- 
age difference is AV = 2 1 14^1. In Fig. [2](a) we used 
AV = 1.0, 5.0 and 10.0 V, generating electric fields of 
0.5, 2.5, and 5.0 V/nm, comparable to experimental val- 
ues. In all experiments^^^^^'^^, the electric displacement 
field in the bilayer is estimated through the dielectric con- 
stant and thickness of the substrate dielectric layer, and 
from the effective potential caused by the carrier doping 
due to the experimental conditions. In this Brief Re- 
port, the electric field was calculated via AV/dbt^ where 
dbt is the distance between the bottom and top gates. 
For the configuration studied here, we verify a decrease 
of the current as AV increases, suggesting the open- 
ing of a band gap that increases with AV. Moreover, 
for fixed gate length and gate voltages, the electric cur- 
rent has a linear dependence with the source-drain volt- 
age, indicating an ohmic contact between the gated and 
non-gated bilayer graphene. This behavior is also ob- 
served experimentally^ From Fig. [2ja) we show 
that we do not reach a zero electric current, even for 
AV = 10.0 V. This result is in agreement with exper- 
imental works that fabricate dual-gate bilayer graphene 
deviceg-^^. 

In Fig. [2Kb) we fixed the Vbg = +2.5 V and Vg = 
—2.5 V, and we present the dependence of the fiowing 
current with the gate length Lgate-, for Lgate = 2 nm, up 
to 10 nm. We also observe a decrease of the current as 
Lgate increases, indicating an increase of the region with 
a gap. Again, the system does not reach a zero current 
even for Lgate = lOnm. 

Since we used open-boundary conditions, with the 
chemical potential fixed by the leads, we can indepen- 
dently vary Vbg and Vtg to control the carriers at the 
scattering region. The presence of the external electric 
field causes a breakdown of the equivalence between the 
two layers. As a result, there is the opening of a band 
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gap and an inter-layer charge transfer. If there were no 
charge transfer, the condition V^g = —Vtg would define 
the charge neutrality point (CNP). However, the pres- 
ence of a Hubbard- like correlation term^^, even at the 
mean-field level, coupled to this charge transfer, leads to 
an overall shift of the bands. As a result, the CNP occurs 
for Vbg —ytg-, as shown in Fig. [3l 




FIG. 3: (Color online) Electronic transport characteristic of a 
bilayer graphene dual-gate FET, for Lgate — ^nm and Vbg = 
-\-2.5V. The current (blue squares) and net charge (black 
circles) at (a) room temperature and at (b) T = 4.5K. 

Thus, we performed calculations where we fixed the 
Vbg at a certain value (+2.5 V) and varied the Vtg. The 
calculated source-drain current and net carrier concentra- 
tion in the bilayer, as a function of the top gate, is shown 
in Figs. [3l^a) and [3jb) for a gate length Lgate = 5 nm, 
for T = 300 and 4.5 respectively. The bias voltage 
was fixed at Vds = 10 mV. The current minimum in 
both curves corresponds to the charge neutrality point. 
For a top gate voltage higher than ^ —3.4 F, electrons 
are being injected in the scattering region, whereas for 
Vtg below this value the holes are being injected. As we 
have no external doping, the chemical potential at the 
gated graphene is at the midgap. Moreover, since the 
conduction and valence bands are very similar, temper- 
ature will create intrinsic carriers but it will not change 
the position of the chemical potential. Thus, there will be 
some charge transfer from the graphene leads (electrons 
or holes) depending on the shift of the mid-gap position, 
at the gated graphene region, relative to the chemical po- 
tential of the leads. Also, in Figs. [3l^a) and[3l^b) we plot 
the carrier concentration per unit area. This presents a 
linear behavior as a function of the variation of the gate 
voltage, with zero net charge at the charge neutrality 
point. This indicates a linear displacement of the bands 
with the gate voltage. There is a notable decrease in the 
current when the temperature is diminished from 300 to 
4.5 i^. This is caused both by variations in the transmit- 
tance with temperature as well as by changes in the tail 
of the Fermi-Dirac distribution [see Eq. (j4])]. For the case 
presented in Fig. [3l the current decreases by more than 
six times. However, the net charge concentration is very 
similar in both temperatures, a result of the invar iance 
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FIG. 4: (Color online) (a) Transmittance for Lgate = 2 nm 
with Vtg = -2.0 y(solid black line), Vtg = -3.1 V (blue 
dashed-dotted line), Vtg = —4.0 V (dashed green line), and 
the nongated case (orange dashed-double-dotted line), (b) 
Transmittance for Lgate = 10 nm with Vtg = —2.0 V (solid 
black line), Vtg = —3.4 V (red dashed-dotted line) and Vtg = 
—4.0 ^/(dashed green line). Inset: Detail of the transmittance 
for Lgate = 10 nm at CNP close to the Fermi energy region. 



of the mid-gap position with temperature, as discussed 
above. 

Figure [H shows the effect of the gate potential in the 
transmittance for (a) Lgate = 2 nm and (b) Lgate = 
10 nm; at 300 K with a fixed back gate voltage of +2.5 V, 
and varying Vtg. Under the same conditions, for larger 
gates the modifications on the transmittance are more 
expressive in comparison to the smaller ones. For each 
Vbg /Vtg configuration the gate potentials causes a low 
transmittance region that is much more pronounced with 
Lgate = 10 nm in comparisson with Lgate = 2 nm. How- 
ever, the transmittance never goes to zero, as can be seen 
on the inset in Fig. lU^b). 

It is important to emphasize that in our calculations 
the electric field is applied only in a finite portion of the 
scattering region through the dual-gate system. As a con- 
sequence, the inter-layer charge transfer caused by this 
electric field occurs only in the gated region. In Fig. [5ja) 
we present the charge induced at CNP, by a dual gate 
with Lgate = 10 nm. In the bottom (top) layer there is 
an excess (lack) of electrons. In Fig. [5fb) we show the 
density of states projected on every four atoms (equiv- 
alent to the unitary cell) along its entire length, from 
the left-hand {z = 0.0 nm) to the right-hand electrodes 
{z = 14.6 nm). Near to the leads (green lines) the pro- 
jected density of states (PDOS) are similar to the PDOS 
of a non-gated pristine bilayer, as expected, whereas di- 
rectly under the gated region (yellow line) the PDOS is 
characteristic of a bilayer with an applied perpendicular 
electric field, where the two peaks located at approxi- 
mately ±12 meV are associated with the presence of a 
"mexican" hat behavior of the bands. However, since the 
gate is spatially finite, even in the presence of an energy 
band gap there is a finite density of states present. This 
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FIG. 5: (Color online) (a) Difference between the non-gated 
and gated local charge density. As can be seen, there is a 
transfer of electrons from the top to the botton layer (units 
of e/bohr^). (b) Projected density of states every four atoms 
(equivalent to the unitary cell) along the z direction. The 
system considered had a Lgate of lOnm and was at the CNP. 



is due to the penetration of the wave function in this for- 
bidden energy region. This results explains the absence 
of an energy gap for transport [e.g.. see Fig. HJb).] 

To understand further the persistent current on the 
device, we look at its behavior as a function of the gate 
length. In Fig. [6l we fix AV = 5 V and Vds = 10 mV, 
and vary Lgate for 300 and 4.5 K. We obtain an ex- 
ponential dependence between the current and the gate 
length, characteristic of a tunneling regime. The tun- 
neling current persists even for the largest gate length, 
Lgate = lOnm. A significant increase of the current 
with temperature can be observed; this is largely caused 
by the broadening of the transport window via the 
[f{E — ijll) — f{E — iir)] term. At the low-bias regime, 
the tail of the Fermi-Dirac function mostly controls the 
current, since there is a small density of states (and hence 
a small number of transport channels) around the Fermi 
level in the electrodes. It is important to remind that in 
our ah initio calculations we retain all the intrinsic prop- 
erties of the device, without the presence of defects, im- 



purity carrier doping, or a substrate dielectric material. 
Therefore, the off-current in our simulations is explained 
by tunneling across the gated region, and strongly de- 
pends on the transport energy window caused by the 
Fermi-Dirac distribution. 

For our nanodevice we estimated the on/off current ra- 
tio, {lon/^off) at room temperature, (300 i^), and low 
temperature, (4.5 K). The on current is calculated for 
the non-gated bilayer graphene with a source-drain volt- 
age of 10 mV and the off currents are all defined at 
the charge neutrality point, with a back gate voltage of 
+2.5 V and the source-drain voltage fixed at 10 mV. 
The largest lon/Ioff happens at Lgate = 10 nm, which 
is around 100 at room temperature, and 1250 at 4.5 K. 

In summary, using ah initio calculations, we study the 
electronic transport properties of a dual-gated bilayer 
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FIG. 6: (Color online) Ids x Lgate curves for AV = 5.0V and 
Vds = lOmV, for 300i^ and 4.5i^. For ah graphs, AV = 2\Vtg\ 
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graphene nanotransistor. Under the action of an external 
electric field, we show that, even for gate lengths up to 
100 A and for room temperature and 4.5 a non-zero 
current is exhibited. The results can be explained by 
the presence of a tunneling regime due the remanescent 
states in the gated region. 
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